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Abstract. The study of sea ice using airborne remote sens- 
ing platforms provides unique capabilities to measure a wide 
variety of sea ice properties. These measurements are use- 
ful for a variety of topics including model evaluation and 
improvement, assessment of satellite retrievals, and incorpo- 
ration into climate data records for analysis of interannual 
variability and long-term trends in sea ice properties. In this 
paper we describe methods for the retrieval of sea ice thick- 
ness, freeboard, and snow depth using data from a multi- 
sensor suite of instruments on NASA’s Operation IceBridge 
airborne campaign. We assess the consistency of the results 
through comparison with independent data sets that demon- 
strate that the IceBridge products are capable of providing a 
reliable record of snow depth and sea ice thickness. We ex- 
plore the impact of inter-campaign instrument changes and 
associated algorithm adaptations as well as the applicabil- 
ity of the adapted algorithms to the ongoing IceBridge mis- 
sion. The uncertainties associated with the retrieval methods 
are determined and placed in the context of their impact on 
the retrieved sea ice thickness. Lastly, we present results for 
the 2009 and 2010 IceBridge campaigns, which are currently 
available in product form via the National Snow and Ice Data 
Center. 


1 Introduction 

The determination of sea ice properties from remote sens- 
ing data has been a long sought-after goal through the uti- 
lization of a wide variety of instruments and field campaigns 
(e.g., Wadhams et al., 1991; Comiso et ah, 2003; Cavalieri et 
ah, 2006; Leuschen et ah, 2008; Haas et ah, 2010). Data on 
sea ice properties from satellite altimetry missions have been 
used to study the behavior of sea ice at regional to global 
scales (e.g., Laxon et ah, 2003; Giles et ah, 2008a; Kwok 
et ah, 2009; Farrell et ah, 2009; Zwally et ah, 2008; Kurtz et 
ah, 2011; Kurtz and Markus, 2012), while data from airborne 
remote sensing missions have been used to validate satellite 
data (e.g., Kurtz et ah, 2008; Connor et ah, 2009; Connor et 
ah, 2012; Laxon et ah, 2013) and identify new ways to extend 
the range of sea ice properties that can be studied using satel- 
lite remote sensing data (e.g., Kurtz et ah, 2009). The sea ice 
properties that remote sensing data sets can retrieve depend 
on the type of instruments used. A key area of interest is the 
use of altimetry data for the study of sea ice thickness. With 
this in mind, a fundamental goal of the ongoing NASA Op- 
eration IceBridge mission is to utilize a suite of instruments 
including radar and laser altimeters to engage in large-scale 
surveys of Arctic and Antarctic sea ice on an annual basis. 
Operation IceBridge is an airborne mission with its primary 
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goal being to create a laser altimetry time series bridging the 
end of the ICESat mission in 2009 and the launch of ICESat- 
2 currently scheduled for 2016. 

The scale of the IceBridge mission in terms of the num- 
ber of instruments flown, areal coverage, and number of ac- 
tive campaigns makes it one of the largest airborne missions 
to date in the polar regions and will enable a new genera- 
tion of research to be conducted using the data sets. A key 
aspect of this study is to describe the retrieval of geophys- 
ical data products from the Level 1 instrument data, which 
can then be used by others wishing to engage in more spe- 
cialized research, as well as to discuss the known limitations 
and uncertainties of the derived geophysical data products. 
Here we focus on three products that can be obtained from 
the IceBridge data set: (1) sea ice freeboard, which is here 
defined as the height of the snow plus sea ice surface above 
the water level to be consistent with the laser altimetry data 
set, (2) snow depth, and (3) sea ice thickness. First, sea ice 
freeboard is the property retrievable by airborne and satellite 
altimeters and provides the capability for sea ice thickness 
studies from the data sets. Moreover, it has been shown to be 
useful as an independent parameterization for energy balance 
studies because of the unique relationship between freeboard, 
sea ice thickness, and the thermal conductivities of snow and 
ice (Kurtz et al., 2011). Second, knowledge of snow depth is 
useful for studies of precipitation trends and variability, melt 
pond coverage, and is also needed for the retrieval of sea ice 
thickness from altimeter data sets. Third, sea ice thickness is 
a long sought-after property which is important for the study 
of a large number of climate variables. Knowledge of this 
property will allow for an assessment of the stability and vari- 
ability of the global sea ice cover, which is a component and 
indicator of global climate change (Lemke et al., 2007). 

Previous studies using IceBridge data over sea ice have 
focused on algorithm development and validation (Farrell et 
al., 2012; Kwok et al., 2012), exploration of instrument be- 
havior and retrieval issues (Kwok et al., 2011), and compar- 
isons to climatological data (Kurtz and Farrell, 2011). How- 
ever, the variable nature of aircraft campaigns means that 
the methods described in previous studies do not necessar- 
ily apply to all of the past and future IceBridge data sets. 
For example, changes to the radar antennas and electronics 
means that algorithms developed for snow depth retrieval for 
the 2009 data set (e.g., Kurtz and Farrell, 201 1) do not auto- 
matically apply to future data sets. Additionally, the distribu- 
tion of leads and collocated laser altimetry data for sea sur- 
face height retrievals along the IceBridge flight tracks also 
changes year to year and means that the uncertainty of the 
retrieved freeboard and ice thickness data is variable along 
the flight paths, rather than a constant value. The addition of 
new instruments and instrument configuration modifications 
allows for retrievals using different and improved method- 
ologies on a year-to-year basis. Essentially, the IceBridge 
data set has and will continue to evolve in time as instru- 
ment changes are made to improve the quality of science re- 



Fig. 1 . Flow chart of processing steps to retrieve sea ice thickness 
from the instrument data. 

trievals from the raw instrument data. In this study, we out- 
line the steps that have been utilized for the retrieval of sea 
ice properties for the 2009 and 2010 data sets, show how 
variations in the retrieval algorithms have changed as a re- 
sult, and discuss the applicability of the new algorithms to 
the 2011 and later data sets. We also present uncertainty esti- 
mates due to variations in the sea ice properties, instrument, 
and inter-campaign algorithm changes. Our goal is to pro- 
vide a baseline reference that describes the main procedures 
for producing a continuous time series of sea ice properties 
and associated uncertainties for the past and future IceBridge 
campaigns. The derived data products described in this work 
can be accessed via the National Snow and Ice Data Center 
at http://nsidc.org/data/idcsi2.html (Kurtz et al., 2012). 

A wide variety of instruments and retrieval methods are 
outlined in this study. Some parameters rely on a single in- 
strument and retrieval method, while others involve multi- 
ple instruments and algorithms. A flow chart describing the 
retrieval of geophysical data products of interest is shown 
in Fig. 1, which can be referenced throughout the text. The 
study is organized as follows: Sect. 2 describes the instru- 
ment data sets. Sections 3, 4, and 5 respectively describe the 
approaches used to retrieve sea ice freeboard, snow depth, 
and sea ice thickness. Each section also includes a discussion 
of error sources due to instrument issues and inter-campaign 
algorithm changes. Section 6 describes the results of the free- 
board, snow depth, and sea ice thickness data sets for the 
2009 and 2010 IceBridge campaigns. The main results of the 
study are discussed and summarized in Sect. 7. 

2 Data sets 

Here we provide an overview of the instruments and data and 
their function in the IceBridge mission. A list of instruments 
and their purpose is provided in Table 1. A map of the 2009 
and 2010 IceBridge flight lines and specific flight numbers 
referred to throughout the text is shown in Fig. 2. 
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Table 1. List of instruments used and their function. 


Instrument 

Geophysical properties retrieved 

Nominal spatial resolution [m] 

Airborne Topographic Mapper (ATM) 

Sea surface height, freeboard 

1 

Snow radar 

Snow depth 

15 

CAMBOT 

Ice type 

0.2 

Digital Mapping System (DMS) 

Ice type 

0.1 

KT-19 infrared pyrometer 

Surface temperature 

130 (for a 1 Hz sampling rate) 




2009 Antarctic IceBridge Sea Ice Flights 

3 Flights from 21 October to 30 October 



Fig. 2. Maps of flight lines for the 2009-2010 Arctic and 2009 Antarctic IceBridge campaigns. Numbers on the map correspond to the 
different flights designated in the text and in Table 5. 


2.1 Laser altimetry data 

The laser altimetry measurements used in this study are from 
the Airborne Topographic Mapper (ATM) (Krabill, 2009a). 
The ATM is a 532 nm wavelength conically scanning laser 
altimeter, combined with a differential GPS system for air- 
craft positioning and an inertial navigation system (INS) 
to measure aircraft orientation. The laser range, GPS posi- 
tion, and INS orientation measurements are used to assign 
three-dimensional geographic coordinates to the point where 
each laser pulse reflects from the surface. The ATM data are 
referenced to the ITRF-2005 reference frame and projected 


onto the WGS-84 ellipsoid. The 15-degree scanner used dur- 
ing the missions yields a measured swath width of approxi- 
mately half of the aircraft’s altitude above the surface. The 
footprint size of each individual elevation measurement is 
~lm, which is set by the laser beam divergence. The sys- 
tem is calibrated using independent ranging measurements 
with the system on the ground, and by overflights of pre- 
surveyed ground areas. Absolute elevation accuracy from the 
ATM is usually about 10cm or better (Krabill et al., 1995) 
with geolocation accuracies of better than 1 m (Schenk et al„ 
1999). Specifically for the IceBridge campaigns, Martin et 
al. (2012) estimate the parameters of the ATM system to be 
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Table 2. Modeled sea surface height data sources. 


Parameter 

Data source 

"geiod 

EGM2008 

"ocean 

TPX06.2 

^load 

TPX06.2 

dearth 

ATM processing (see text) 

"pressure (AP) 

ECMWF Interim reanalysis 


(1) 74 cm horizontal accuracy, (2) 6.6 cm vertical accuracy, 
and (3) 3 cm vertical precision. 

Using ground test calibration data, Kwok et al. (2012) 
identified a bias in the ATM data for low amplitude signals 
due to the leading edge fit procedure used for elevation re- 
trievals. The bias will impact the freeboard retrievals when 
points with a low signal strength are used. We have applied 
the polynomial fit correction factor described in Kwok et 
al. (2012) to correct for these biases. No correction was ap- 
plied for points with a signal strength higher than 950 (the 
ATM signal strength is not radiometrically calibrated and is 
thus reported as dimensionless data numbers) since the poly- 
nomial fit diverges beyond this point. A different laser sys- 
tem was operated during the 2010 Arctic and 2009 Antarctic 
campaigns, thus we have applied a similar correction gen- 
erated from ground calibration data specific for these cam- 
paigns. The correction, h t _ c orr (in unit of meters), is an 8th 
order polynomial fit of calibration data provided by the ATM 
team that is added to the ATM elevation data: 

/t e - corr = 1-356 x 10“ 26 r s 8 - 1.51483 x 10“ 22 r s 7 (1) 

+7.48991 x 10“ 19 r s 6 - 2.16621 x 10“ 15 r s 5 
+3.97857 x 10“ 12 r s 4 -4.61175 x 10“ 9 r s 3 
+3.17998 x 10“ 6 r 2 - 0.001 18755 r s + 0.2, 

where r s is the ATM received signal strength. This polyno- 
mial was chosen to correct for biases to the level of points 
with a received signal strength of 1 100, which was the modal 
signal strength during the 2010 Arctic campaign. A constant 
correction of 0.008 m is added to points with a received sig- 
nal strength greater than 2500 to correspond to the calibra- 
tion data set and account for divergence of the polynomial 
fit beyond this level. The range calibration data and derived 
h e-con correction factor are shown in Fig. 3a, and histograms 
of the received signal strength values are shown in Fig. 3b. 
The signal strength has a mean of 93 1 and standard deviation 
of 387 for the 2009 campaign, and a mean of 1064 and stan- 
dard deviation of 397 for the 2010 campaign. For the 2009 
campaign, the mean value of the bias correction over sea ice 
covered points is 0.005and 0.023 m over leads. For the 2010 
campaign, the mean value of the bias correction is 0.003 m 
over sea ice covered points and 0.032 m over leads. 



Received signal strength 


Received signal strength 


Fig. 3. (a) Dependence of received signal strength on measured 
range to a fixed target for the 2010 ATM instrument. Black dots 
are the mean elevation bias for each unit of signal strength. The 
black line is the subsequent correction factor that was added to the 
ATM elevations used in this study for a particular signal strength, 
(b) Histogram of ATM received signal strength values for the 2009 
and 2010 Arctic campaigns. 


2.2 Snow radar 

Snow depth was measured using a frequency-modulated, 
continuous-wave (FMCW) radar that utilizes dechirp on re- 
ceive (Leuschen, 2009; Panzer et al., 2013). Mixing a de- 
layed, attenuated receive signal with a copy of the transmit 
signal returns a beat frequency equal to the product of the 
chirp rate (bandwidth divided by pulse length) and the two- 
way delay time to the target. Dining the 2009 campaign, the 
radar transmitted 2.5 to 7 GHz with a 270 ps pulse length and 
2 kHz pulse repetition frequency. During the 2010 campaign, 
the radar transmitted 2.0 to 6.5 GHz with a 250 ps pulse 
length and 2 kHz pulse repetition frequency. Four hardware 
presums (also referred to as coherent averaging or stack- 
ing) are performed prior to writing data to disk with an 
additional four presums done in post-processing. Processed 
radar data are provided in 33 second segments along the 
flight path, spanning 4.2 km for the nominal aircraft speed 
(460 km hr-'). 

The snow radar system measures the return radar signal as 
a function of time, which is scattered from the area illumi- 
nated beneath the aircraft. Snow depth is determined by de- 
tecting the snow-air and snow-ice interfaces within the radar 
waveform and multiplying the time separation between the 
interfaces by the speed of light within the snow pack. Details 
of the snow depth retrieval process are provided in Sect. 4. 
Given the 4.5 GHz bandwidth, the theoretical free-space ver- 
tical range resolution is approximately 5 cm. For detection of 
the snow layer over sea ice, the vertical range resolution im- 
proves to approximately 3 cm due to speed of light variations 
in the snow pack. However, during post-processing, the data 
are Hanning windowed to reduce range sidelobes, thus de- 
grading the vertical range resolution by a factor of 2. At the 
nominal flight altitude (460 m) the radar has a footprint size 
of 1 1 m across track (dictated by the pulse-limited footprint 
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size) and 14.5 m along track (dictated by the synthetic aper- 
ture formed). The footprint size reduces to the diameter of the 
first Fresnel zone of ~ 8 m when coherent (specular) returns 
dominate the receive signal. For specular returns to occur, the 
Rayleigh criterion must be satisfied over the first two Fresnel 
zones of the radar active area (Carsey, 1992). The Rayleigh 
criterion is defined as surface roughness < |, where X is the 
radar wavelength. For the snow radar, the Rayleigh criterion 
is satisfied for surface roughness values less than ~ 1 cm. 
Such smooth ice is expected to be only rarely encountered 
over the snow-covered sea ice portions of the survey areas 
and thus diffuse rather than specular reflections dominate the 
returns by the radar. Therefore, we take the synthetic aper- 
ture dictated footprint size to be the size of the radar active 
area for our snow depth retrievals. 

2.3 Visible imagery 

We use two sets of aerial photography data for identifying 
morphological features on the sea ice in conjunction with the 
radar and laser altimetry data. The use of the aerial photog- 
raphy data was found to be essential to distinguish surface 
types, which is needed for our sea ice freeboard retrieval 
process. This also required knowledge of the geolocation of 
each pixel within the aerial photography data sets. The Digi- 
tal Mapping System (DMS) (Dominguez, 2009) is available 
for all flight campaigns after the first Arctic campaign in 
2009. The DMS is a digital camera that takes high resolu- 
tion (~ 10 cm) natural color (RGB) photographs that span 
the length of each IceBridge flight line. The angular resolu- 
tion of the system is sufficient to cover the full scan width 
of the laser altimeter measurements. The data are geolocated 
and the latitude and longitude of each image pixel were de- 
termined using metadata parameters provided in each image 
file. 

For the 2009 IceBridge Arctic campaign, the DMS instru- 
ment was not available, however, the CAMBOT (Continu- 
ous Airborne Mapping By Optical Translator) imaging sys- 
tem was available (Krabill, 2009b). The CAMBOT images 
are taken once every five seconds along each flight path us- 
ing a Canon EOS Digital Rebel XTi camera. For the 2009 
campaign, the CAMBOT images were not overlapping and 
there is ~ 300 m of track that is not covered between each 
image. The images are time tagged and the geolocation of 
the center point, height of the aircraft above the surface, 
and aircraft heading are provided from the aircraft trajec- 
tory information based on the time tag. The angular dimen- 
sions of the CAMBOT images for the 2009 Arctic campaign 
are 58.12° x 40.71°. The aircraft altitude and angular dimen- 
sions were used to determine the pixel size for each CAM- 
BOT image in the 2009 data set. The pixel size was combined 
with the image center point location, aircraft heading, and a 
standard coordinate rotation to geolocate the pixels of each 
CAMBOT image. However, the time-tagging procedure for 
the CAMBOT data set was typically found to be valid only 


Cambot Image 0013, March 31 si, 2009 



Fig. 4. Example of a CAMBOT image of sea ice with ATM ele- 
vation (h corr) measurements (colored circles) manually overlain so 
that features and elevations correspond. 


to within ±1 second, which introduces geolocation errors for 
each image. To refine the geolocation of each image, we have 
designed software to manually align each CAMBOT image. 
The software overlays the more accurately geolocated ATM 
elevation data onto each image, and then allows the operator 
to adjust the CAMBOT geolocation by manually changing 
the geolocation of the center point of each image. This was 
done until topographic features such as ridges and leads coin- 
cided in both the images and the ATM data. Figure 4 shows 
an example image of a manually aligned CAMBOT image 
with ATM data overlain. 


3 Sea ice freeboard retrievals 

Various methods for the retrieval of freeboard from the ATM 
airborne laser altimetry data have been described in the lit- 
erature (e.g., Kurtz et al„ 2008; Connor et ah, 2009, 2012; 
Kwok et ah, 2012). Here, we describe a new approach utiliz- 
ing geolocated aerial photography and a lead discrimination 
algorithm to maximize the quality and number of laser al- 
timeter data points used to determine the sea surface height. 
We have found this to be an important component of the free- 
board retrieval step due to the loss of a large percentage of the 
laser altimeter returns over thin ice-covered and open-water 
leads (subsequent IceBridge campaigns have utilized a sec- 
ond 3-degree off-nadir ATM laser system, which minimized 
the loss of laser returns from leads). This loss of data oc- 
curs due to specular reflection of the laser pulse away from 
the receiver when insufficient surface roughness elements are 
present to cause diffuse scattering. There is still a sufficient 
number of returns available over leads to determine sea sur- 
face height, but this loss introduces two problems: (1) the 
low (and variable) number of returns over leads compared 
to the large number of returns from the surrounding sea ice 
make the returns from leads difficult to uniquely distinguish 
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in surface elevation histograms. (2) The loss of laser altime- 
try data over thin sea ice contributes to an undersampling of 
these regions and thus a biased freeboard result when taking 
an areal average. To resolve these issues for more accurate 
retrieval of sea ice freeboard we use a combination of pho- 
tography and laser altimetry as described in this section. 

The primary product given by the ATM laser altimeter is 
the surface elevation referenced to the WGS84 ellipsoid, h e . 
However, the more geophysically useful parameter we wish 
to retrieve from this product is the height of the snow plus ice 
surface above sea level, termed the freeboard, fb. The conver- 
sion of elevation data into sea ice freeboard is accomplished 
by subtracting out the instantaneous sea surface height h ss h 
from each elevation measurement: 

fb = h e /? ssh - 

The instantaneous sea surface height at a given point in space 
and time can be written as (Chelton et al., 2001) 

/( ssh — ^geoid 4~ /( tides 4~ “pressure 4" /(dynamic, (2) 

where /7 ge oid is the geoid, /(tides is the contribution of tidal 
forces, /? pressure is the effect of atmospheric pressure loading, 
and /(dynamic is the dynamic topography of the ocean surface. 
We first estimate the sea surface height along the flight tracks 
using modeled estimates for the h ge 0 id, /(tides, and /(pressure 
terms. Data sources used in the initial estimation of the sea 
surface height terms are shown in Table 2. The geoid model 
is the Earth Gravitational Model 2008 (EGM2008) (Pavlis et 
al., 2008). The /(tides component can be further decomposed 
into 4 terms: 

/hides — /(ocean 4“ /( load + /(earth 4“ /(pole, 

where /(ocean is the ocean tide, /i load is the load tide, /(earth 
is the solid earth tide, and /( po ie is the pole tide. The ATM 
elevation product is provided in the ITRF standard tide-free 
system from which the solid earth tide has been removed. A 
standard latitude dependent correction is also applied to the 
data used here to place the data in a mean tide system relative 
to the WGS84 ellipsoid. The pole tide is a small amplitude 
(<~ 2 cm), long wavelength tidal component caused by os- 
cillations in the Earth’s rotation axis (Wahr, 1985), it has not 
been included here since its impact on freeboard determi- 
nation is negligible (less than 1 mm of elevation difference 
over a length scale of 100 km). The TPX06.2 tide model is 
used to estimate the /(ocean and /(load components (Egbert and 
Erofeeva, 2002). The isostatic response of the sea surface to 
atmospheric pressure loading, /(pressure , can be rewritten as 

_ A P 

/(pressure = , 

Pv/g 

where A P is the difference between the surface air pressure 
at the local point from the instantaneous mean surface air 
pressure over the ocean, p w = 1024kgm -3 the density of sea 


water, and g — 9.8 ms -2 the gravitational acceleration (Chel- 
ton et al., 2001). In this data set, A P is calculated as the dif- 
ference between linearly interpolated local surface pressure 
data provided by the European Center for Medium-Range 
Weather Forecasts (ECMWF) ERA-Interim data set and an 
average global mean surface air pressure of 101 300 Pa. 

As a first step, we remove variations in the instantaneous 
sea surface height along the flight tracks through the subtrac- 
tion of modeled estimates for the /i ge oid, /(tides, and /(pressure, 
resulting in a corrected elevation /( corr : 

/(corr — he (/(tides 4“ /(geoid T /(pressure)- 

The modeled sea surface height parameters do not at this 
point have sufficient accuracy for the useful retrieval of free- 
board given that a 1 cm freeboard error corresponds to a 
~ 10 cm error in sea ice thickness. The periodic identifica- 
tion of local sea surface elevations for use as tie points is 
required to achieve the desired centimeter level accuracy for 
freeboard retrievals. After the removal of the modeled sea 
surface height parameters, the resulting freeboard is calcu- 
lated as 

fb = /(corr — Zssh, (3) 

where z ss [, is the locally determined sea surface elevation 
with respect to the /z corr elevation data set. The determina- 
tion of z ss h at each point along the flight track follows from 
the retrieval of individual sea surface height observations de- 
scribed in detail below. 

A set of sea surface elevation estimates, h t p , are first found 
through extraction of the /( corr elevation data identified over 
leads in the IceBridge visible imagery. For all but the Arctic 
2009 data set, an automated lead detection algorithm called 
Sea Ice Lead Detection Algorithm using Minimal Signal 
(SILDAMS) (Onana et al., 2012) that utilizes the geolocated 
IceBridge DMS imagery is used. SILDAMS applies a min- 
imal signal transformation on DMS pixel brightness values 
to carry out a localization around low pixel intensities that 
correspond to leads. The transformed outputs are within a 
uniform dynamic variability over the set of numbers from [0, 
1] for a variety of input image pixel intensities. This allows 
specified thresholding of the transformed outputs to be set 
for the identification and classification of three different lead 
classes corresponding to open water areas, grease ice/nilas, 
and newly frozen leads with non-snow-covered grey ice (see 
Fig. 5). The thresholds are set through manual classification 
of select lead images within a flight segment and the SIL- 
DAMS thresholds are adjusted to optimally match the manu- 
ally classified ice types. The SILDAMS output for each DMS 
image is visually examined to ensure the quality of the out- 
put, the thresholds are maintained along the flight line until 
visual examination reveals an error in the output, which oc- 
curs at discrete time periods due to camera adjustments to 
account for changing lighting conditions. Each of these lead 
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Elevation 



Fig. 5. Example of a DMS image of sea ice with ATM elevation (/icon ) measurements (colored circles) overlain. The right-hand image shows 
results from the SILDAMS algorithm that is used to distinguish between (1) open water leads (red), (2) nilas and grease ice (green) and (3) 
non-snow-covered grey-white ice (blue). 


Table 3. Summary of radar system parameters impacting the receive signal power and signal-to-noise ratio from the snow-air interface for 
the 2009 and 2010 campaigns. 


Parameter 

Symbol in text 

2009 

2010 

Transmit power 

Pt 

5 mW 

20 mW 

Gain 

G 

lOdBi 

8 dBi 

Center wavelength 

X 

0.063 m 

0.071 m 

Compressed pulse length 

r 

2.22 x 10 -10 s -1 

2.22 x 10“ 10 s“ 1 

Time bandwidth-product 

Tbp 

60.85 dB 

60.51 dB 

Number of integrations 

N, 

16 

16 

System temperature 

T 

290 K 

290 K 

System noise figure 

F 

21.32 dB 

19.08 dB 


classes has an ice thickness of less than 30 cm (World Meteo- 
rological Organization, 1970) and should thus have an eleva- 
tion that falls within 3 cm of the local mean sea level making 
them suitable for use as sea surface tie points. To account for 
biases in the use of thin ice areas as sea surface tie points, 
we subtracted 0.005 and 0.02 m from the sea surface height 
estimates derived from the ATM elevation data over grease 
ice/nilas and grey ice, respectively. These values were chosen 
to correspond to expected freeboards of snow-free ice types 
defined in the World Meteorological Organization nomencla- 
ture. 

For the 2009 Arctic campaign, leads and ice type were 
categorized through visual inspection of CAMBOT images 
(Krabill, 2009b). Visual lead identification, rather than an au- 


tomated approach using SILDAMS, was also used due to the 
complexity of the CAMBOT images caused by failure of the 
mechanical camera shutter and uneven exposure in some por- 
tions of the images. By matching the location of the leads in 
the visible imagery with those of the ATM h corT elevations, 
the local sea surface elevation and freeboard of each ATM 
elevation point was retrieved in the same manner as leads 
identified with SILDAMS. Bias removal due to the use of 
thin ice types as sea surface tie points was carried out in the 
same manner as described for the SILDAMS approach. We 
note that even if leads are present in an area, the number of 
sea surface elevation estimates may be much more limited in 
number than those from an equivalent area of returns over 
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sea ice due to the loss of ATM data from smooth specularly 
reflecting surfaces (as can be seen in Figs. 4 and 5). 

The spatial size of each individual ATM laser footprint is 
1 m for the nominal flight altitude (460 m) of the IceBridge 
data set. To account for geolocation uncertainties (0.74 m), 
as well as the presence of mixed lead/sea ice data within the 
ATM returns, the requirement of a minimum 1 m lead buffer 
was set for each ATM return over a lead for it to be used as a 
sea surface estimate. This is necessary because mixed lead- 
sea ice returns are prevalent due to the highly backscattering 
sea ice portion of the surface; as such, laser returns near the 
edge of a lead are not representative of the actual sea surface 
elevation but rather that of the sea ice within the laser foot- 
print. Thus, the set of sea surface elevation estimates, /z tp , are 
taken from the h con elevation data set where leads are iden- 
tified, and the lead is found to extend an additional 1 m in all 
directions beyond the edge of the laser footprint. 

Within a given area, the mean sea surface elevation, /z tp , is 
derived from the set of sea surface elevation estimates, /z tp , 
that were determined using the combined ATM and visible 
imagery lead detection method described previously. /? tp is 
calculated using a procedure described in detail below. First, 
all values of /z tp within ±250 m from the center point are first 
combined into a histogram with a 2 cm bin size to initially 
determine the local sea surface height. A value of ±250 m 
has been chosen to span the width of an individual DMS 
or CAMBOT image. We expect the local sea surface height 
to be constant over a length scale corresponding to the first 
mode baroclinic Rossby radius of deformation, which is on 
the order of 10 km for latitudes greater than 60° and less than 
-60° (Chelton et al., 1998). The Rossby radius is associated 
with the length scale at which oceanic eddies form; these 
eddies can cause local inhomogeneities in the sea surface 
height. The length scale of ±250 m used for the initial sea 
surface height determination is much less than the Rossby 
radius, and thus we expect the actual sea surface height to 
be constant over this chosen length scale. The standard de- 
viation of the /? tp estimates in this local area is then a mea- 
sure of the uncertainty of the individual sea surface height 
estimates due to instrument noise. An analysis of histograms 
of ATM elevations over known flat surfaces, including sep- 
arate cases with open water and flat snow-covered sea ice 
showed that the elevation distributions are Gaussian in shape 
with a standard deviation of ~ 10 cm or less. This has also 
been observed in separate studies by Farrell et al. (2012) 
who reported a standard deviation of 5 cm over a flat, snow- 
covered refrozen sea ice lead, while Connor et al. (2012) 
reported a standard deviation of 10 cm over lead areas. We 
thus ideally expect the distribution of all h tp points within 
the length scale less than the Rossby radius to be similar in 
shape and width. However, a variety of error sources ranging 
from laser mounting biases, geolocation errors, misidentifi- 
cation of lead returns in the visible imagery, and errors in the 
/z ssh data sources lead to deviations from this ideal scenario. 
These errors also preclude the use of the mean value of the 


local /z tp observations and standard error analysis techniques 
from being used to determine /z tp and its associated error. The 
misidentification of lead returns within the combined ATM 
and photography data causes the largest impact, which, when 
present, causes the presence of secondary and higher modes 
in the histogram of /z tp points. In determining /z tp , we wish to 
use points corresponding to the mode with the lowest eleva- 
tion. This is accomplished through the use of the centroid of 
a Gaussian fit function to the histogram of h tp . The follow- 
ing conditions have been imposed to determine whether the 
fitted Gaussian function is of high enough quality for use in 
determining h tp : 

1 . erfit < 0. 1 1 m, 

2. x 2 <0.015, 

3. N >40, 

where ofr is the standard deviation of the Gaussian fit, / 2 is 
the reduced chi-square goodness-of-fit, and N is the number 
of h t p points used to construct the histogram for the Gaus- 
sian fit. In cases where a multi-modal distribution of /z tp is 
observed, the above parameters will not be satisfied for the 
initial Gaussian fit. If this occurs, an iteration is then per- 
formed by discarding the largest h tp elevation point in the set, 
performing another Gaussian fit, and retesting the fit param- 
eters. The iteration is repeated until the conditions for the fit 
parameters are satisfied, at which point /z tp is determined. If 
the conditions for the fit parameters are not met, then the sea 
surface height and freeboard are not calculated. See Fig. 6a 
for an example of a case where the initial Gaussian fit param- 
eters were satisfied on the first iteration, and Fig. 6b for an 
example case of a multi-modal distribution where the itera- 
tion produced a fit of the first mode only. 

This procedure produces a set of discrete sea surface 
height observations (or tie points) along the flight track. 
Fig. 7 maps the locations of these sea surface height observa- 
tions for the Arctic campaigns. The availability of these sea 
surface height observations is non-uniform along each flight 
track due to ice conditions and dynamics. In particular, few 
leads were found to occur in the compact multi-year ice ar- 
eas of the Arctic. This creates uncertainties in the retrieval of 
freeboard since the accuracy of the sea surface height obser- 
vations decreases with distance due to remaining uncertain- 
ties in the removal of the modeled sea surface height param- 
eters shown in Eq. 2. Therefore, to construct the along-track 
freeboard profiles and determine the uncertainty of each free- 
board estimate (which is inherently variable along the flight 
path due to uneven lead spacing), the sea surface height along 
the full flight track was interpolated using an ordinary krig- 
ing approach (Cressie, 1993). Although the basin-scale sea 
surface height field is non-isotropic and non-homogeneous 
(due to variability in the components of Eq. 2), we assume 
that the covariance of the sea surface height field is a ho- 
mogeneous function of distance. The interpolation is then 
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a) b) 




Fig. 6 . Histogram of ATM elevations (/i C orr) over lead areas identi- 
fied using the SILDAMS algorithm. The Gaussian fits used to de- 
termine the instantaneous sea surface height are also shown for two 
cases: (a) a Gaussian fit which satisfied the fit criteria on the first 
iteration, and (b) a Gaussian fit which satisfied the fit criteria after 
the iteration procedure was completed. 

performed by minimizing the error variance of the interpo- 
lated column vector of sea surface heights for all tie points 
within ±200 km of each observation point along the flight 
path. Briefly, z S sh is calculated at a given point as 

z ssh = W T z„, (4) 


where d xy is the distance between points x and y, L is the 
correlation length of the sea surface height observations (z 0 ) 
along the flight segment (determined from the location where 
the autocorrelation sequence falls below e 1 for each flight), 
a z is the standard deviation of the sea surface height observa- 
tions for the flight, and e ss h is the uncertainty of the local sea 
surface height observations. From the previously described 
method for determining /? tp containing at least 40 observa- 
tions and a standard deviation of 0. 1 1 m, the uncertainty for 
the sea surface height observations, e ss h, should theoretically 
be better than 0.1 1 /a/40 = 0.017 m. However, additional er- 
rors such as unresolved aircraft pitch and roll errors will 
likely make this number higher since each sea surface height 
elevation data point, h tp , is not statistically independent. We 
assume that each of the h tp observations are statistically inde- 
pendent and thus take e SS h = 0.058/(^" =1 exp(-jj))- 5 me- 
ters (with a maximum value of e ss h = 0.058 m) based on ob- 
servations from the 2009 Arctic data set where 0.058 m is the 
standard deviation of the differences of all h tp measurements 
within ±5 km (the Rossby radius where we expect the sea 
surface height to be constant) of each individually observed 
lead. Using a length scale smaller than ±5 km did not change 
the maximum value of e ss h = 0.058 m suggesting that this is 
the minimum uncertainty of a single sea surface height ob- 
servation due to instrumental limitations. Lastly, the sea ice 
freeboard is then determined by combining Eqs. 3 and 4. The 
uncertainty in sea surface height and therefore freeboard at 
each point is given by 


where z S sh is the interpolated sea surface height at a specific 
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where C is the covariance matrix and C oy is the covariance 
between the interpolation point y and the observation point 
o, and /x is a Lagrange multiplier. The elements of the co- 
variance matrix are assumed to be a homogeneous function 
of distance, which is modeled as a Gaussian process with the 
respective elements taken to be 


Cxy — e ss h ' 




Qualitatively, the local sea surface height is determined by 
weighting the sea surface height observations as a function of 
distance, with nearby h tp points receiving the largest weight. 
The uncertainty decreases as the number of sea surface ob- 
servations increases and the distance from each sea surface 
observation decreases, approaching er ss h = m in an area 
with multiple sea surface height observations. The uncer- 
tainty approaches that of the background sea surface height 
field standard deviation, a z , as the distance from the nearest 
sea surface height observation becomes large. 

This method is used to calculate the freeboard and its asso- 
ciated uncertainty at each individual ATM data point (~lm 
spatial resolution ) where the local sea surface height is de- 
termined. However, in order to combine the freeboard data 
with the snow depth estimates (described in the next section) 
for the retrieval of sea ice thickness, we place the high res- 
olution freeboard estimates onto the same spatial resolution 
as the snow depth measurements. This is done by averag- 
ing the individual freeboard data to a 40 m resolution profile 
along the flight track centered at the snow radar footprints. 
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However, due to the loss of ATM data over surfaces with low 
backscatter (leads and non-snow covered thin ice, such as 
those examples illustrated in Figs. 4 and 5), an adjustment to 
the profile data is required to account for the undersampling 
of the freeboard of thin ice. This adjustment is performed 
by filling ATM data gaps using the SILDAMS output as a 
proxy freeboard data set. The freeboard of points correspond- 
ing to open water, / ow , were set to zero. The mean freeboard 
of points corresponding to “thin ice” in the SILDAMS al- 
gorithm were set to be / t ; = 0.005 m corresponding to the 
WMO freeboard (and ice thickness) for grease ice/nilas. Sim- 
ilarly, points corresponding to the “grey ice” classification in 
the SILDAMS algorithm were set to be / g ; = 0.02 m. Based 
on these estimates, we determine an adjusted average free- 
board through the combined use of the ATM and SILDAMS 
freeboard outputs over a 40 m area. The ATM and SILDAMS 
data were placed on a 1 m resolution grid. The 40 m resolu- 
tion adjusted freeboard, fb a dj- was found using the following 
relation: 

(^ow/ow A ZZ ti ,/ti + ^gi/gi T n jeefb) 

fbadj = 7 7 , 

(now + T zZgi -j- H{ce) 

where n x corresponds to the number of open water, thin ice, 
grey ice, and sea ice grid points in the 40 m area detected by 
SILDAMS, and fb is the mean freeboard from the ATM data 
within the area. This adjusted freeboard represents the profile 
of the mean freeboard along the IceBridge flight line where 
data gaps over low backscattering ice and water surfaces have 
been taken into account. 

The methodology applied for the retrieval of freeboard is 
independent of location and thus is used to retrieve freeboard 
for both the Arctic and Antarctic IceBridge campaigns. How- 
ever, for the 2009 Arctic campaign no DMS data were avail- 
able for the automated identification of leads. CAMBOT data 
are available for all of the 2009 Arctic campaign with the 
exception of ~ 3/4 of the April 5th flight where mechanical 
problems prevented photographs from being taken. Thus, it 
was possible to retrieve freeboard for the 2009 Arctic cam- 
paign, but estimates could only be processed to the fb level, 
rather than to fb a dj • 

4 Snow depth retrievals 

Methods for the retrieval of snow depth from FMCW radar 
operation over sea ice have been described in previous stud- 
ies (e.g., Galin et al., 2011; Kwok et al., 2011). Retrieval 
methods for the IceBridge snow radar have been described 
by Kurtz and Farrell (2011), and Farrell et al. (2012). The 
method of Kurtz and Farrell (201 1) is used to retrieve snow 
depth for the 2009 IceBridge campaign in this study. This 
method also forms the basis of a new method described in 
this study for the 2010 and future Arctic IceBridge cam- 
paigns. We note that the snow depth retrieval method de- 
scribed in this document may be applicable only for the Arc- 


tic region. For Antarctic sea ice, observations have shown 
that due to the complicated snow morphology that is often 
found (Massom et al., 2001), the strongest return in a radar 
signal may be reflected from an undetermined point within 
the snow layer rather than the snow-ice interface (Giles et 
al., 2008b; Willatt et al., 2010) . The presence of flooded sea 
ice will also complicate the retrieval of snow depth from the 
radar system. Therefore, we consider the retrieval of snow 
depth from the snow radar system for the Antarctic region to 
be a research area that requires in-depth study to understand 
the results and determine the associated uncertainties. 

Briefly, the retrieval algorithms for the snow radar system 
detect the snow-air and snow-ice interfaces within the radar 
waveform and determine the snow depth by multiplying the 
time separation between the interfaces by the speed of light 
within the snow pack. Comparison of the radar signal with 
the in situ observations of Farrell et al. (2012) have shown 
that the strongest peak in the return is expected to corre- 
spond to the snow-ice interface, but that the snow-air inter- 
face does not often correspond to a distinct peak as may be 
expected from theoretical arguments. Factors such as instru- 
ment noise and response, volume scattering, snow density 
variations, and surface roughness features all combine to cre- 
ate a complex signal that deviates from the ideal theoretical 
signal where two distinct peaks are expected to correspond to 
the snow-air and snow-ice interfaces. Therefore, the method 
used here to detect the air-snow interface is empirical, but 
has been shown to produce accurate results in comparison to 
in situ data (Farrell et al., 2012). 

Modifications to the radar hardware and installation be- 
tween IceBridge campaigns preclude the use of the snow 
depth retrieval algorithm for the 2009 campaign from be- 
ing applied to the 2010 and future IceBridge campaigns. For 
the 2009 campaign, the sub-optimal transmit/receive antenna 
isolation required reducing the transmit power and attenuat- 
ing the receiver front end to avoid saturation. Both measures 
resulted in decreased loop sensitivity. This was evident in 
the 2009 data as the detected snow-air interfaces were near 
the noise floor, and the snow depth retrieval algorithm of 
Kurtz and Farrell (2011) was specifically designed to work 
under these conditions. In 2010, the radar was rebuilt and 
the operating frequencies were changed from 2.5-7 to 2.0- 
6.5 GHz. By improving the transmit/receive antenna isola- 
tion (approximately 35 dB in 2009 to better than 45 dB in 
2010), the transmit power was increased from 5 mW in 2009 
to 20 mW in 2010 and the front-end receiver attenuation was 
reduced thus improving the receiver noise figure and reduc- 
ing the thermal noise level. Table 3 shows the differences 
between the radar system parameters for the 2009 and 2010 
campaigns. The radar equation for a pulse-limited altimeter 
(Raney, 1998) outlines the basic set of parameters that influ- 
ence the received power from the snow-air interface: 

PtG 2 X 2 O 0 7TCXTBP 

Pr = 5 , ( 6 ) 

(47 tR) 3 
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Fig. 7. Map of lead locations identified in the aerial photography data sets for (a) 2009 and (b) 2010. Black boxes represent areas where 
leads were identified in the aerial photography and the red points represent leads which contained sufficient ATM measurements suitable for 
the retrieval of the instantaneous sea surface height. 


Table 4. Summary of all uncertainty components propagated in the calculation of the final IceBridge sea ice thickness uncertainty. 


Parameter 

Symbol in text 

Uncertainty 

Source 

Snow density 

°Ps 

100 % 

Climatological variability 

Sea ice density 

ffp, 

10 if 

Estimated variability from in situ measurements 

Snow depth 

°h s 

0.057 m 

Estimate from in situ comparison 

Freeboard 

f 

Variable (~ 1 cm to 30 cm) 

Kriging error (see Sect. 3) 


Table 5. Mean freeboard, snow depth, and sea ice thickness for the 
2009-2010 IceBridge flights. The flight numbers correspond to the 
flight lines mapped in Fig. 2. The mean freeboard and thickness 
values are reported only using data where the freeboard uncertainty 
is < 0. 1 m, the mean ice thickness uncertainty is in parentheses. The 
mean snow depth is reported for the entire flight track where snow 
depth is > 0.05 m and less than the derived freeboard. 


where P r is the received power, P, is the transmit power, G 
is the antenna gain, k the wavelength, ct° the backscattering 
coefficient, r is the compressed pulse length (equal to the re- 
ciprocal of the radar bandwidth), 7bp is the range time band- 
width product, and R is the range to the target. The signal- 
to-noise ratio (SNR) is given by 


SNR = 


P,G 2 X 2 0 O 7TctT bp N i 

(4ttR) 3 IcTBF 


(7) 


Flight 

date 

Flight 

number 

Mean 

freeboard [m] 

Mean snow 
depth [m] 

Mean sea ice 
thickness [m] 

3/31/2009 

i 

0.42 

- 

- 

4/2/2009 

2 

0.48 

0.31 

2.50 (0.74) 

4/5/2009 

3 

0.45 

0.21 

2.54 (0.65) 

4/21/2009 

4 

0.67 

0.38 

3.82 (0.75) 

4/25/2009 

5 

0.51 

0.28 

3.71 (0.66) 

10/30/2009 

6 

0.59 

- 

- 

3/23/2010 

7 

0.52 

0.34 

3.39 (0.78) 

3/26/2010 

8 

0.32 

0.26 

2.91 (0.67) 

4/2/2010 

9 

0.18 

0.15 

1.77 (0.64) 

4/5/2010 

10 

0.51 

0.24 

3.64 (0.64) 

4/12/2010 

11 

0.53 

0.29 

3.62 (0.68) 

4/19/2010 

12 

0.38 

0.23 

2.58 (0.58) 

4/20/2010 

13 

0.48 

0.26 

3.22 (0.72) 

4/21/2010 

14 

0.43 

0.25 

3.00 (0.59) 


where Nj is the number of coherent integrations, k is the 
Boltzmann constant, T is the system temperature, B is the 
bandwidth, and F is the system noise figure. To demonstrate 
the improved SNR between the two campaigns, we use Eq. 7 
and take a ratio for the two campaigns (assuming a constant 
backscattering coefficient, flight altitude, center wavelength, 
and the values shown in Table 3) that yields an increased 
SNR of ~ 4.9 dB for the received signal from the snow-air 
interface for the 2010 campaign. The ~ 4.9 dB improvement 
in SNR can be primarily attributed to the increase in transmit 
power. 

The improved radar design and subsequent increased 
signal-to-noise ratio for the 2010 data and future data sets 
means that modifications to the snow depth retrieval algo- 
rithm of Kurtz and Farrell (2011) are required. To account 
for these changes in the radar configuration, in the 2010 
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campaign we base the detection thresholds for the air-snow 
interface on the relative radar return power. We use the 2009 
radar data set as the basis for the algorithm changes since 
the in situ study of Farrell et al. (2012) demonstrated the 
successful detection of the snow-air interface by the algo- 
rithm, and thus provides us with a reference data set of the 
return power values from the snow-air and snow-ice inter- 
faces. The location of the snow-ice interface is still defined 
as the largest maxima in the radar signal below the air- 
snow interface. The return power in dB is calculated here as 
P — 10 x log 10 ( 77 ^^ where /n data is the magnitude of the 
echogram data and /V cc h 0 is the number of points in an indi- 
vidual radar waveform. To improve the SNR and reduce the 
possibility of radar speckle from being misidentified as an in- 
terface layer, the return power data were first incoherently av- 
eraged (incoherent averaging rather than coherent averaging 
is done since phase information was not originally provided 
in the 2009 data set) over 40 points corresponding to a length 
scale of ~ 40 m . The return power values for the snow-air 
interfaces identified in the 2009 data set were found to vary 
between P™" = — 4.0dB and P s n l ; “ = —0.5 dB, with a mean 
value of P s '™“' = — 2.35dB, and these same conditions are 
used to locate the snow-air interfaces in the 2010 data set. 
However, the data from the snow radar are presently not ra- 
diometrically calibrated and thus vary based on the parame- 
ters of the instrumentation, the flight parameters (Eq. 6 ), and 
the processing methodologies used to create the final data 
values. Thus, the return power values for the 2010 IceBridge 
campaign are not directly equivalent to those from the 2009 
data set. 

Given the differences between the radar system in 2009 
and 2010 , differences due to flight operation considerations 
(such as altitude), and lack of radiometric calibration for the 
return power values, a calibration adjustment for each indi- 
vidual radar echogram is therefore necessary to determine 
the final threshold power range where we expect the snow- 
air interface to be located. First, we define scale, m , and off- 
set, b, parameters, which are determined for each echogram 
by simultaneously solving 


yi — mx i + b , ( 8 ) 


y 2 = mx 2 + b , (9) 

where y\ = 2.25 dB, which is the average power of the 
snow-ice interface observed in the 2009 Arctic data set, x\ is 
the average power of the snow-ice interface for the desired 
radar echogram, y 2 — —5.0 dB is the mean noise level, which 
is defined as the average power of the first 100 bins located 
at least 5 m above the point of maximum power (assumed to 
be in the air above the snow pack) in the 2009 Arctic data 
set, and X 2 is the mean noise level, which is similarly de- 
fined as the average power of the first 100 bins located at 


least 5 m above the point of maximum power for the desired 
radar echogram. Here we define the maximum power of the 
adjusted snow-air interface as P™* = mP™* + b, similarly 
the minimum power for the adjusted snow-air interface is de- 
fined as P™" = mP™ + b. We use these conditions to de- 
termine the region of the radar return where we expect the 
snow-air interface to be located. The search for the snow-air 
interface thus occurs when the return power level, P s - a , for 
the snow-air interface satisfies the following requirements: 


-> pinin 
— r s—a ’ 


( 10 ) 


> pmin 
— r s—a ’ 


(ii) 


where P s - a is the mean power in the six range bins that 
follow the first point with power > P™". These conditions 
identify the beginning portion of the return from the snow-air 
interface by finding the point where the minimum expected 
return power from the interface is reached (first condition) 
and when volume scattering from within the snow pack is 
also present (second condition). Once these conditions are 
satisfied the snow-air interface is defined following the ob- 
servations described in Farrell et al. (2012) to be either the 
location of the first local maxima in the return or the point 
of continuous power increase. The snow-air interface is thus 
defined here as the first point where one of the following con- 
ditions is satisfied: 


1 . The point of continuous power increase of the radar re- 
turn from the snow-air interface is found. The point 
of continuous power increase is defined here as the 
point where the radar return power begins to contin- 
uously increase (i.e., -jf- > 0 ) until the maximum ex- 
pected snow-air interface power, P s “ a * is reached. 

2. If a local maxima occurs greater than m P s ™ ean + b and 
a above the adjacent points (to eliminate random noise 
from being misidentified as a maxima), where a is the 
standard deviation of the radar noise level, then this 
point is chosen as the snow-air interface location. 

Examples of two radar waveforms from the 2010 data set il- 
lustrating cases where a distinct local maxima is not found 
(condition 1) and is found (condition 2) are shown in Fig. 8 . 
The figure shows that the theoretical case of an easily dis- 
tinguishable maxima corresponding to a reflection from the 
snow-air interface is not observed in the data. Thus, the 
power range between P™ “ and P s ™“ is used to identify the 
distinct region in the return waveform where the snow-air 
interface is located. 

An iteration of the above method is performed when the 
difference between the echogram mean noise level and the 
echogram mean snow-ice interface power is greater than 
twice the mean difference observed in the 2009 campaign 
data set (equal to 14.5 dB). The iteration is performed by 
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Fig. 8. (a) Example of calculated snow-air (red line) and snow-ice (black line) interfaces for a radar file from the 2010 data set. (b) Example 
radar waveform illustrating the range of power values where the snow-air interface is located. In this case, the snow-air interface corresponds 
to the first point of continuously increasing power rather than a distinct peak. The snow-ice interface is the point of maximum power, (c) 
Example radar waveform illustrating the case where a distinct peak is found for the snow-air interface. 


rescaling the threshold power values: y\ is set equal to the 
mean echogram power for the initially determined snow-air 
interface and v ’2 = P ™™ . The initial scaling provides an ap- 
proximation for the parameters m and b by scaling the return 
power values by the snow-ice interface power and the system 
noise, however as previously demonstrated, changes to the 
snow radar system have reduced the noise level in successive 
campaigns so the initial scaling is only an approximation. 
The iteration scales the return power values to the snow-air 
interface allowing for more accurate estimates of m and b to 
be made when there is a large difference between the mean 
noise level and snow-ice interface power. 

The snow depth is then found by calculating the delay time 
between the snow-air and snow-ice interfaces and multiply- 
ing this difference by the speed of light within the snow pack, 
Csnow- The speed of light in the snow pack was taken to follow 
the relation between dry snow density and dielectric constant 
given by Tiuri et al. (1984) as 

ed — 1 + 2p s , 

c 


where c is the speed of light in vacuum and p s is the snow 
density (in ^|j). The snow density was taken to be 320 
following the climatological mean snow density provided in 
Warren et al. (1999). Once the snow-air and snow-ice inter- 
faces are identified, we apply a locally weighted robust linear 
regression at a 40 m length scale to reduce the impact of out- 
liers in the final determination of the snow-air and snow-ice 
interface locations. This effectively reduces the spatial reso- 
lution of the retrieved snow depths to 40 m. 

As described in Kurtz and Farrell (2011), the behav- 
ior of the radar over leads and the apex of steep pressure 
ridges must also be accounted for in the snow depth retrieval 
method. Lead areas were flagged using the method described 
in Kurtz and Farrell, (2011) and the snow depth is set to zero 
to correspond to the negligible snow cover on open water 
and newly frozen leads. We also discard data where the sig- 
nal is too low for the retrieval of snow depth. Based on the 
analysis of coincident ATM laser data used to identify lo- 
cations where there is large disagreement between the laser 
and radar snow-air interfaces, the conditions for discarding 
the data due to insufficient signal strength were set as 

Ps-i > — 1.5m+Z?[dB], (12) 
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Fig. 9. Example radar waveform over a ridge. The point of maxi- 
mum power in the radar waveform is within the expected range for 
the snow-air interface and does not reach the minimum required 
threshold, P™?. 


Ps-i > -1.5m + fc[dB], (13) 

where P s _,- is the mean power in the three range bins that 
follow the estimated snow-ice interface, and m and b are 
found by solving Eqs. 8 and 9. An example radar waveform 
is shown in Fig. 9 over a sea ice ridge that was identified in 
the ATM and DMS data sets. This loss of signal typically 
occurs over the apex of steep pressure ridges, however, com- 
parison with in situ data (Farrell et al., 2012) showed this 
data loss did not limit the ability of the snow radar to cap- 
ture the mean snow depth of a region since the snow sur- 
rounding the ridges was still retrieved. Additionally, Kurtz 
and Farrell (2011) demonstrated the difference between the 
snow radar for the full 2009 data set is within 0.3 cm of the 
climatology of Warren et al. (1999) and provides further sup- 
port that the mean snow depth over a variety of ice types 
is captured. If the snow depth data set from the snow radar 
was not representative of the mean snow depth over the full 
range of ice types it could potentially lead to errors in the 
retrieval of sea ice thickness. How representative the snow 
radar data is on basin scale will continue to be examined as 
future studies with coincident in situ survey lines conducted 
in 2011 and 2012 are completed. In total, 16% of the snow 
depth data was discarded for the 2009 campaign, and 21 % 
was discarded for the 2010 campaign. 

In addition to discarding data due to insufficient signal 
strength, it is also necessary to filter data that have been neg- 
atively impacted by other physical constraints. As shown by 
Kurtz and Farrell (2011) it is necessary to discard data col- 


lected with warm surface temperatures due to the presence 
of liquid water, which changes the dielectric properties of the 
snow pack. Similar to the study of Kurtz and Farrell (2011), 
we discard data when the surface temperature is greater than 
-5 °C, this is the temperature threshold identified in observa- 
tions by Barber et al. (1995) where large dielectric changes 
in the snow pack begin to take place. When available, the 
surface temperature is determined from the IceBridge KT-19 
infrared radiation pyrometer data set (Shetter et al., 2010). 
When instrumental observations are unavailable, we use the 
thermodynamic sea ice model of Kurtz et al. (2011) forced 
with ECMWF meteorological data to determine the surface 
temperature. Lastly, flight altitude variations also need to be 
taken into account. The nominal IceBridge flight altitude 
is 460 m, but is variable due to mission requirements and 
weather considerations. A maximum altitude of 540 m for 
the 2009 campaign and 675 m for the 2010 campaign is set 
for the retrieval of snow depth. These altitude limits corre- 
spond to the point where the expected snow-air interface re- 
turn power is near to the system noise level for the respective 
campaigns given a l/R 3 return power dependence. 

The method applied for the retrieval of snow depth dif- 
fers between the 2009 and 2010 seasons. The method ap- 
plied to the 2009 data set was validated through comparison 
to in situ data, which gives us a measure of confidence in 
the snow depths for that season (Farrell et al., 2012), as well 
as a demonstrated agreement with the climatology of War- 
ren et al. (1999) (Kurtz and Farrell, 2011). However, no in 
situ surveys were conducted in 2010, which precludes a thor- 
ough assessment of the 2010 data set. Multiple overflights 
of in situ surveys were conducted in subsequent years (as 
of this writing, in 2011 and 2012), which will allow more 
detailed analysis of the retrieval methods and associated un- 
certainties described in this study. To assess the quality of 
the snow depth retrievals from the 2010 data set we compare 
the snow-air interface retrieved from the radar to that derived 
from the ATM surface elevation data. Figure 10 shows an ex- 
ample of a radar echogram from the 2010 data set with the 
chosen snow-air (red line) and snow-ice interfaces (black 
line). Also shown is the ATM corrected elevation (/7 corr ) data 
(green line), which should ideally correspond to the snow- 
air interface detected by the radar. However, the radar data 
are referenced as a relative distance from the snow radar an- 
tenna, rather than to a standard geodetic coordinate system 
such as the WGS84 ellipsoid used for the ATM data. There- 
fore, we take a mean difference between the ATM elevation 
(averaged to a 40 m resolution) and the chosen snow-air in- 
terface to construct the plot. The figures show the ability of 
the tracking algorithm to identify the snow-air interface, and 
also show the behavior over ridges and thin ice. We used the 
coincident ATM data to assess the ability of the algorithm 
to determine the quality of the snow-air interface. A proba- 
bility distribution showing the correlation between the ATM 
data and the identified snow-air interface for each snow radar 
echogram is presented in Fig. 11. The figure illustrates that 
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Fig. 10. 2010 IceBridge ATM and snow radar comparison for (a) level flying by the aircraft and (b) non-level flying corresponding to aircraft 
pitch and roll changes of 3 degrees beginning near the 1 km along-track distance portion of the flight. The green line is the ATM snow-air 
surface, the red line is the calculated snow-air interface from the radar data, the black line is the calculated snow-ice interface. The correlation 
between the ATM and snow radar picked snow-air interface is 0.75 for the level flying case and 0.19 for the non-level flying case. 


the correlation is quite variable. The maximum correlation of 
the ATM and radar data is near 1 with a mode around 0.55 
demonstrating the ability of the snow-air interface algorithm 
to track the snow surface in certain cases. The modal correla- 
tion of 0.55 is also similar to the results of Farrell et al. (2012) 
who observed a maximum correlation of 0.7 when compar- 
ing the IceBridge derived snow depths with in situ data. 

The observation that many echograms have only moder- 
ate to low correlations can be explained due to a combina- 
tion of uncertainties in the snow depth retrievals (such as 
due to the finite range resolution of ~ 5 cm), geolocation 
uncertainties in the snow radar data, and aircraft motion. 
Low variability of the snow-air interface location can also 
reduce the correlations as the ATM and radar system noise 
will dominate in these cases. Additional uncertainties in the 
radar data comparison arise due to the uncorrected lever arm 
between the GPS antenna, to which the INS data is refer- 
enced to, and the monostatic phase center of the radar. This 
influences the comparison in two ways: (1) it leads to er- 
rors in the vertical range compensation, which are used to 


account for changes in the aircraft altitude and (2) it intro- 
duces a 5-20 m offset between the geolocations of the snow 
radar and ATM. While we expect that averaging to a 40 m 
length scale reduces the problem of geolocation errors of this 
magnitude for sea ice thickness retrievals using the radar and 
ATM freeboard data, it nonetheless is a source of error if 
one seeks to carry out a one-to-one matching of the ATM 
data over the snow radar footprint such as done by Kwok 
et al. (2011). Aircraft pitch and roll introduces errors in the 
aircraft altitude correction that results in residual features in 
the vertical dimension of the radar data when they are com- 
pared to the ATM elevation data that have been corrected for 
aircraft pitch and roll. Figure 10 shows two example radar 
echograms where good correlation between the ATM surface 
elevation and radar retrieved snow-air interface over small 
length scales can be seen visually, but Fig. 10b demonstrates 
the impact of uncorrected aircraft motion on the snow radar 
data, which can be clearly seen when compared to the motion 
corrected ATM elevation data. The pulse-limited processing 
and large beamwidth of the antenna (45° by 45°) will still 
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Fig. 11. Probability density function of the correlation between the 
ATM-derived snow-air interface and the picked radar snow-air in- 
terface over each ~ 4 km radar echogram for the 2010 Arctic cam- 
paign. 


produce a return from nadir even when the aircraft is not 
level (typical flight operating parameters do not have roll an- 
gles greater than 20 degrees), however the uncorrected lever 
arm between the GPS and radar antennas introduces errors 
in the aircraft altitude correction, which appear as vertical 
variations (Fig. 10b). Since we define the snow depth as the 
relative difference between the snow-air and snow-ice inter- 
faces, and the pitch and roll typically does not change sig- 
nificantly over the 40 m retrieval length scale, we expect that 
aircraft motion will have little impact on our snow depth re- 
trievals. However, these pitch and roll fluctuations introduce 
a randomly varying, but continuous source of error in the cor- 
relation between the ATM and snow radar snow-air interface 
locations. For the 2009 Arctic campaign, the standard devia- 
tion of aircraft pitch was 0.25 degrees and the standard devia- 
tion of roll was 2.6 degrees. Given the magnitude of the pitch 
and roll variations across the flight lines, we use the corre- 
lations with the ATM freeboard data as qualitative evidence 
that our snow-air tracking algorithm is capable of accurate 
retrievals under ideal conditions (as evident in the right-hand 
tail of the correlation histogram shown in Fig. 1 1). However, 
the combined ATM-snow radar data set cannot currently be 
taken as a quantitative assessment of the snow-air interface 
retrieval quality. Future planned comparisons with in situ 
data sets are needed to determine a quantitative estimate of 
the snow depth retrieval quality. 


5 Sea ice thickness retrievals 
5.1 Thickness determination 

Sea ice thickness, hi, is calculated using the corresponding 
40 m scale freeboard and snow depth data sets as input for 
the hydrostatic balance equation: 


1 Pw ,, Pw Ps , .. 

hi = — fbadj — h s , (14) 

Pw Pi Pw Pi 

where fb a( jj is the freeboard (described in Sect. 3), h s is 
the snow depth (described in Sect. 4), p w is the density of 
sea water, p ; - is the density of sea ice, and p s is the density of 
snow. p w and p,- are taken to be 1024 -^4 and 915 ^4, which 
are derived from the result of numerous field measurements 
summarized by Wadhams et al. (1992). p s is taken to be 320 
following the climatological values compiled by Warren 
et al. (1999). 

5.2 Error analysis 

The error in the sea ice thickness retrieval (excluding the neg- 
ligible contribution of errors due to variations in sea water 
density) can be written as 
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where aj tf , Oh s , cr Ps , and a Pi are the uncertainties of 
the ice thickness, freeboard, snow depth, and densities of 
snow and ice, respectively. A summary of the uncertainty 
terms is provided in Table 4. Density uncertainties are taken 
from previous in situ measurements described in the litera- 
ture: a Ps is estimated to be 100 kg m -3 based on the variabil- 
ity of p s in the climatology of Warren et al. (1999), a Pi is es- 
timated here to be 10 kg m -3 , which represents the expected 
range of densities for sea ice between 0.3 and 3 m thick (Ko- 
vacs, 1996) as well as the typical range of sea ice densities 
from observations (Wadhams et al., 1992). We note that the 
density of sea ice may be dependent on ice type, with the den- 
sity of multi-year ice typically being lower than that of first 
year ice due to the presence of air bubbles (Wadhams et al., 
1992; Alexandrov et al., 2010). The mean density of sea ice 
has also been found to be variable with respect to the amount 
of ice that is above the water line, as well as the sea ice salin- 
ity and temperature (Timco and Frederking, 1995; Eicken et 
al., 1995). The use of a density value that is based on a spe- 
cific ice type or geographical region will be investigated in 
future studies. The uncertainty in the freeboard retrieval, Oh f , 
is variable along the flight path as described in Sect. 3. The 
uncertainty in snow depth is not well constrained at the time 
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of this writing. Uncertainty in the snow depth occurs due to a 
variety of factors including the finite range resolution of the 
radar, density uncertainties, and uncertainty in the detection 
of the snow-air and snow-ice interfaces. Considering only 
uncertainties in range resolution and speed of light varia- 
tions, Kwok et al. (201 1) calculated an error of 3.5 to 5 cm for 
snow depths between 10 and 70 cm. This represents the min- 
imum expected uncertainty in snow depth derived from the 
radar. Following the discussion in Sect. 4 on issues involving 
detection of the snow-air interface, we expect that uncertain- 
ties in interface detection will play a much more significant 
role in the uncertainty of the snow depth retrieval. The com- 
bined uncertainties due to speed of light variations and inter- 
face detection uncertainties mean that a robust error estimate 
for an individual snow depth measurement should be derived 
by computing the standard deviation of differences between 
coincident radar and in situ data sets. However, spatial off- 
sets between the IceBridge data and in situ measurements de- 
scribed in Farrell et al. (2012) preclude a robust error assess- 
ment from being done, and no other coincident in situ survey 
data exist for the 2009-2010 time period. For the purposes of 
this study, we estimate the uncertainty using the difference 
between the IceBridge data and in situ snow depth described 
in the study by Farrell et al. (2012). The mean difference 
between the survey and IceBridge data set was 0.8 cm, and 
there were 50 independent observations (40 m spatial reso- 
lution for each radar observation over a ~ 2 km survey line). 
Assuming the IceBridge and in situ data are unbiased, the 
eiror is calculated as ct/, s = 0.8 x V50 = 5.7 cm. While this 
is a simplistic error analysis, it is a useful estimate of the 
expected uncertainty since it incorporates the important ef- 
fect of interface detection issues. However, we await further 
refinement of this uncertainty value through the analysis of 
multiple coincident in situ data surveys conducted during the 
2011 and 2012 campaigns. Improvements to the snow radar 
and increased signal-to-noise ratio over the 2009 to present 
time period may also lead to campaign specific uncertainty 
estimates. 

The uncertainty in ice thickness is calculated using Eq. 16, 
using the uncertainties for each variable described above and 
in Table 4. The sea ice thickness uncertainty is thus a vari- 
able quantity, in particular due to the variable uncertainty in 
the freeboard retrievals. It is also important to note that the 
uncertainties in the freeboard are correlated over long length 
scales due to the procedures for determining the sea surface 
height, therefore, spatial averaging of the data to reduce the 
uncertainty will need to take this into account. 

We provide an estimate of sea ice thickness over a 40 m 
length scale to provide the highest resolution available from 
the data. However, since sea ice is a rigid body the hydro- 
static balance equation does not necessarily apply for each 
of the ~40m resolution freeboard and snow depth mea- 
surements taken along the flight path. An analysis of Arc- 
tic sea ice by Forsstrom et al. (2011) found an uncertainty 
of 0.5 m associated with non-local hydrostatic balance. Sim- 


ilarly, Doble et al. (201 1) showed that the ice thickness error 
in the use of the hydrostatic balance equation was 0.4 m for 
averaging length scales of 300 m. The dependence of this er- 
ror on length scale was observed to be dependent on the ice 
type, the error in level ice decreased faster with increasing 
averaging length than the error of deformed ice. The use of a 
40 m length scale used here for calculating the sea ice thick- 
ness will thus have an additional error term due to local non- 
hydrostatic balance effects. This error is important to con- 
sider when comparing the IceBridge data to small-scale mea- 
surements, but will become minimal when averaging over 
length scales greater than the size of an individual ice floe 
since the large-scale ice pack is in hydrostatic equilibrium. 

6 Results for 2009 and 2010 

In this section we show the results of the retrieved freeboard, 
snow depth, and sea ice thickness data for the 2009-2010 
IceBridge campaigns. For the purposes of this section we 
only use data where the uncertainty in the retrieved free- 
board is less than 0.1 m. Maps of the sea ice thickness and 
snow depth from the Arctic campaigns are shown in Figs. 12 
and 13, while mean values are shown in Table 5. Histograms 
of the sea ice thickness and snow depth distributions for FY 
(first year) and MY (multi-year) ice are shown in Fig. 14. 
The maps show the varying spatial coverage of the data sets 
with increased coverage during the 2010 campaign due to a 
higher number of flight surveys. The snow depth maps show 
nearly full coverage over the flight tracks while the ice thick- 
ness maps show considerably less coverage due to (1) the 
absence of visible imagery for the detection of leads during 
much of Flight 3; (2) the absence of snow depth retrievals 
for Flight 1 (due to the higher transmit power and associated 
receiver saturation issues for this flight); and (3) the sporadic 
locations of leads necessary for accurate (a/, ; < 0.1 m) free- 
board retrievals. The increased number of flights since 2009 
has led to the availability of more sea ice thickness data in 
each campaign. 

The maps of snow depth show spatial gradients that are 
largely consistent with the expected pattern of snow distribu- 
tion from observations (e.g., Warren et al., 1999). The deep- 
est snow cover is found north of Greenland and the Canadian 
Archipelago with lower snow depths towards the Beaufort 
Sea region. As described in Kurtz and Farrell (2011), first 
year ice regions were found to have significantly less snow 
than multi-year ice regions, which is not reflected in the cli- 
matology of Warren et al. (1999) since the region was largely 
covered by multi-year ice when the observations were taken. 
This can also be seen in the snow thickness distributions 
in Fig. 14. For the 2009 campaign, the mean radar derived 
snow depth over the first year ice regions (16.4 cm) was ob- 
served to be 15 cm (48 %) less than the mean snow depth 
over multi-year ice regions (31.6 cm), while for the 2010 
campaign the mean snow depth over first year ice (13.2 cm) 
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Fig. 12. Maps illustrating the along-track snow depth derived for the Arctic IceBridge campaigns for (a) 2009 and (b) 2010. The red line 
corresponds to the multi-year ice boundary from the AMSR-E mask. The numbers correspond to the flight lines in Table 5. 


la) 


8 




2 3 4 5 

Sea ice Thickness 

Fig. 13. Maps illustrating sea ice thickness for the Arctic IceBridge campaigns for (a) 2009 and (b) 2010. The red line corresponds to the 
multi-year ice boundary from the AMSR-E mask. 



was 14 cm (52 %) less than the mean snow depth over multi- 
year ice (27.4 cm). For multi-year ice, the mean snow depths 
are close to the climatology of Warren et al. (1999) (Kurtz 
and Farrell, 2011). The mean snow depth over first year ice 
is also broadly consistent with the mean snow depth from 
the AMSR-E snow depth on sea ice data product (Cava- 
lieri et al., 2004). Gridding the IceBridge snow radar data 
to the same 12.5 km polar stereographic grid as AMSR-E re- 
sults in mean differences (radar— AMSR-E) of 1.3 cm for the 
2009 campaign and -1.2 cm for the 2010 campaign. Farrell 
et al. (2012) found a mean difference of 0.8 cm between the 


IceBridge derived snow depths and a 2 km in situ line of data 
during the 2009 campaign. We note that these comparisons 
of the radar derived mean snow depth values are not meant to 
supplant a more detailed study of the data sets and associated 
errors, but to show the snow depth estimates presented here 
are consistent with other independent data. 

The ice thickness maps also show large-scale spatial gra- 
dients that are consistent with the expected pattern from re- 
cent observations (e.g., Haas et al., 2010; Kwok et al., 2009; 
Laxon et al., 2003; Wadhams et al., 2011; Laxon et al., 
2013) and model results (e.g., Schweiger et al., 2011). The 
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Fig. 14. Histograms of the ice and snow thickness distributions for 
the 2009 and 2010 campaigns separated into FY and MY ice com- 
ponents. The FY ice thickness data for the 2009 campaign are pri- 
marily front the Nares Strait region and are not necessarily repre- 
sentative of the FY ice thickness distribution from the Arctic basin. 


thickest ice is observed north of Greenland and the Canadian 
Archipelago where the convergence of the sea ice pack due to 
the mean circulation pattern of the Arctic Ocean leads to the 
prevalence of heavily deformed and ridged thick ice. Away 
from the coastal areas and towards the Beaufort Sea region 
lower ice thickness values are observed. The thickness of first 
year ice has been observed to be significantly lower than that 
of multi-year ice, and this is reflected in the IceBridge data 
set. For the 2010 campaign (very few thickness observations 
are available for first year ice in the 2009 data set), the mean 
thickness of the first year ice area is 1.99 m while the mean 
thickness of the multi-year ice area is 3.36 m. 

We note that establishing the absolute accuracy of the Ice- 
Bridge data set remains an ongoing effort, but here we com- 
pare the consistency of the results so as to provide support 
for the differences in inter-campaign retrieval methods as de- 
scribed in this document. In order to compare the consistency 
of the 2009 and 2010 results we have first gridded the Ice- 
Bridge data to a 25 km polar stereographic grid. There were 
38 grid points with spatial overlap corresponding mainly to 
Flights 2 and 3 for 2009, and the repeat Flights 12 and 14 for 
2010. The mean thickness for the 2009 campaign overlap- 
ping grid points was 2.65 m, while the mean thickness for the 
2010 campaign was 2.61 m. While we do expect differences 
due to interannual variability as well as ice motion, the con- 


sistency of the results is encouraging. Flights 4 and 10 further 
demonstrate the inter-campaign consistency of the retrieved 
results. These two flights were not exact repeats and thus not 
present in the previous comparison, but instead were closely 
spaced flight tracks that flew underneath the orbit of the En- 
visat satellite. Table 5 shows that we observed a large dif- 
ference of 0.14 m in the mean snow depth between the 2009 
and 2010 flights, but we also observed a large corresponding 
change of 0.16 m in the measured mean freeboard. The total 
difference in mean ice thickness was 0.18 m. This difference 
is small considering the large mean differences in freeboard 
and snow depth and the independence of the two data sets. 

Lastly, as a measure of the consistency within a campaign 
we compare overlapping 25 km grid points within the 2010 
campaign data set (the 2009 campaign had only three over- 
lapping grid points, which is not sufficient for a comparison). 
There were 58 overlapping grid points with snow depth mea- 
surements, the mean difference was 0.024 m with a standard 
deviation of 0.066 m. This compares well with the expected 
uncertainty of 0.057 m reported earlier and expected differ- 
ences (such as ice motion or snow fall events) due to the ~ 1 
month time period over which the campaign was conducted. 
For ice thickness, there were 13 overlapping grid points with 
a mean difference of 0.22 m and a standard deviation of 
0.78 m. This also compares well with the mean uncertainty of 
the grid points, which was derived to be 0.70 m, and also ex- 
pected differences due to time offsets. Taken together, these 
results show that there is consistency in the IceBridge snow 
depth and sea ice thickness results in comparison to indepen- 
dent data sets as well as inter-campaign and intra-campaign 
consistency. 


7 Conclusions 

Here we have presented a description of the methods used to 
retrieve sea ice freeboard, thickness, and snow depth using 
a synthesis of instrument data from the IceBridge airborne 
platform. The derived products for the 2009 and 2010 Ice- 
Bridge campaigns are currently available for public usage at 
http://nsidc.org/data/idcsi2.html. The methods described in 
this manuscript will also be applied for the retrieval of sea 
ice properties from the 2011 and 2012 campaigns and will 
also be made available through the National Snow and Ice 
Data Center. The inclusion of geolocated aerial photography, 
laser altimetry data, and snow radar data are all expected to 
be available throughout the IceBridge mission and thus the 
retrieval methods developed here will serve as a core pro- 
cessing methodology for future IceBridge campaigns. The 
inclusion of additional instruments is expected to lead to im- 
proved retrievals in subsequent years. For example, the KT- 
19 infrared pyrometer was included in the 2012 IceBridge 
campaign and was set to use a sampling rate of 0.3 Hz, which 
was higher than the 1 Hz sampling rate used in previous cam- 
paigns, this has the potential to improve the retrieval of sea 
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ice freeboard by providing a method to discriminate sea ice 
leads when the sun is near or below the horizon and thus 
aerial photography is not available. This is an area of future 
research and any new retrieval method using the KT-19 data 
set will first be compared to the standard freeboard retrieval 
method that uses laser altimetry data and aerial photography 
described in this document. 

Overall, we have shown there is consistency within the 
IceBridge snow depth and sea ice thickness results, and also 
shown there is consistency in comparison to other indepen- 
dent data sets. Linking the IceBridge data set to the sea ice 
thickness record, in particular to the records from satellite al- 
timetry data such as ICESat, CryoSat-2, and ICESat-2, will 
require detailed knowledge of uncertainties in the data set. 
The evaluation and refinement of the accuracy of the Ice- 
Bridge data set remains an ongoing effort especially with re- 
gard to the determination of the uncertainty at the 40 m length 
scale of the retrievals as well as the determination of any 
residual biases in the basin-scale data. This will be studied 
in future comparisons to other instruments and coincident in 
situ data from subsequent IceBridge campaigns. 

Sea ice thickness and snow depth data derived using the 
core methodology described in this document have the po- 
tential to support a wide variety of research activities. A ma- 
jor ongoing effort is the connection of the IceBridge data set 
with the record of sea ice thickness retrievals from ICESat 
and CryoSat-2 to ensure that a consistent record of Arctic 
sea ice thickness is achievable from the satellite altimetry 
record. The results have also been used to support seasonal 
forecasting of the sea ice minimum and analyze the impact 
that initialized ice thicknesses have on the predictive capa- 
bilities of models (Lindsay et al., 2012). There is also po- 
tential for the evaluation and improvement of current snow 
depth retrieval methods such as those currently used from 
passive microwave data (Markus and Cavalieri, 1998; Cav- 
alieri et al., 2012). Many other research efforts will benefit 
from the availability of such a data set and it is our goal to 
ensure that a consistent, and high quality, set of data prod- 
ucts are produced to support the ongoing research goals of 
the cryospheric community. 
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